A System of Differential Equations for the Airy Process 



Craig A. Tracy 

Department of Mathematics 
University of California, Davis, CA 95616 
email address: tracy@math.ucdavis.edu 

Harold Widom 

Department of Mathematics 
University of California, Santa Cruz, CA 95064 
email address: widom@math.ucsc.edu 

Abstract 

The Airy process r — ► A T is characterized by its finite-dimensional distribu- 
tion functions 

Pr(^ n <d,...,A Tm <£m). 

For m = 1 it is known that Pr (A T < £) is expressible in terms of a solution 
to Painleve II. We show that each finite-dimensional distribution function 
is expressible in terms of a solution to a system of differential equations. 

I. Introduction 

The Airy process r — > A T , introduced by Prahofer and Spohn [6], is 
the limiting stationary process for a certain 1 + 1-dimensional local random 
growth model called the polynuclear growth model (PNG). It is conjectured 
that the Airy process is, in fact, the limiting process for a wide class of 
random growth models. (This class is called the 1 + 1-dimensional KPZ 
universality class in the physics literature [5].) The PNG model is closely 
related to the length of the longest increasing subsequence in a random 
permutation [2]. This fact together with the result of Baik, Deift and Jo- 
hansson [3] on the limiting distribution of the length of the longest increasing 
subsequence in a random permutation shows that the distribution function 
Pr (A T < £) equals the limiting distribution function, i<2(£)> 01 the largest 
eigenvalue in the Gaussian Unitary Ensemble [7]. i*2 is expressible either as 
a Fredholm determinant of a certain trace-class operator (the Airy kernel) 
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or in terms of a solution to a nonlinear differential equation (Painleve II). 
The finite-dimensional distribution functions 

Pr(A T1 <^,...,A Tm <U) 

are expressible as a Fredholm determinant of a trace-class operator (the 
extended Airy kernel) [4, 6]. It is natural to conjecture [4, 6] that these 
distribution functions are also expressible in terms of a solution to a system 
of differential equations. It is this last conjecture which we prove. For 
m = 2 this conjecture was also proved, in a different form, by Adler and van 
Moerbeke [1]. 



II. Statement 



The Airy process is characterized by the probabilities 

Pr (A T1 < . . . , A Tm < U) = det (J - K), 
where K is the operator with m x m matrix kernel having entries 

Kij(x,y) = L ij (x,y)X(t j , 00 )(y) 



and 



Lij(x,y) = < 



e -z (Ti-Tj) _|_ z ) Ai(y + z) dz if i > j, 



r ■ ' - T J> Ai(x + z) Ai(y + z)dz if i < j. 



We assume throughout that r\ < ■ ■ ■ < r m , and think of -ftT as acting on the 
m- fold direct sum of L 2 (a, oo) where a < min £j. 

To state the result we let R = K (I — K)^ 1 and let A{x) denote the 
mxm diagonal matrix diag (Ai(x)) and X(x) the diagonal matrix diag (Xj(x)), 
where Xj = X(g. j00 ). Then we define the matrix functions Q(x) and Q(x) by 

Q = (I-K)~ 1 A, Q = AX(I-K)- 1 

(where for Q the operators act on the right). These and R(x, y) are functions 
of the £j as well as x and y. We define the matrix functions q, q and r of 
the ij only by 



Qij = QiMi) 



We always interpret Rij(x, £,■) as the limit Rij(x, These quantities are inde- 

pendent of our choice of a. 
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Finally we let r denote the diagonal matrix diag (jj ) . 

Our differential operator is V = where dj = d/d^j, and the 

system of equations is 



Here the brackets denote commutator, £ denotes the diagonal matrix diag 
and 6 the matrix with all entries equal to one. 

This can be interpreted as a system of ordinary differential equations if 
we replace the variables £1, £ m by & + £,..., £ m + £, where £1, . . . , £ m 
are fixed and £ variable. Then V = d/d£, and the £j are regarded as 
parameters. 

To get a representation for det (7 — K) observe that 



where the last factor denotes multiplication by the diagonal matrix with all 
entries zero except for the j th , which equals 5(x — We deduce that 

dj log det (J — K) = -Tr (/ - K)' 1 djK = Rjj^j, 

Hence V log det(7 — K) = Tr r, and so it follows from (3) that 

V 2 log det (7 — K) = — Tr q@q 

since the trace of [r, r] equals zero. This gives the representation 



Here the determinant is evaluated at . . . , £ m ) and in the integral £ + ?? 
is shorthand for (£i + 77, . . . , £ m + 77). 

If m = 1 the commutators drop out, q = q, equations (1) and (2) are 
Painleve II and these are the previously known results. 




Sq + 2qeqq-2[T, r]q, 
q£ + 2qq@q-2q[T, r], 
-q&q+ [r, r]. 



(1) 
(2) 
(3) 



Vr 



dj K = —L8j, 



(4) 
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III. Proof 



The proof will follow along the lines of the derivation in [7] for the case 
m = 1. There the kernel was "integrable" in the sense that its commutator 
with M, the operator of multiplication by x, was of finite rank. The same 
was then true of the resolvent kernel, which was useful. But now our kernel 
is not integrable, so there will necessarily be some differences. 

With D = d/dx we compute that 
[D, K\ij = -Ai(x) Ai(y) Xj(y) + L tj (x, 5(y - + (n - Tj ) K tj {x, y). 
Equivalently, 

[D, K] = -A(x) A(y) X(y) +L5 + [r,K], 

where S = J2j &j-> multiplication by the matrix diag (5(x — and L is the 
operator with kernel Lij(x, y). (For clarity we sometimes write the kernel of 
an operator in place of the operator itself.) To obtain [D, R] we replace K 
by K — I in the commutators and left- and right-multiply by p = (I — K)^ 1 . 
The result is 

[D, R] = -Q(x)GQ(y) + RSp+ [r, p}. 2 (5) 
We have already defined the matrix functions Q and Q and we define 

P = (J - K)- l A\ u = (Q, Ai) = J Q(x) Ai(x) dx. 

It follows from (5) and the fact that r and A commute that 

Q' = P-Q@u + RSQ+[t, Q}:> (6) 

Next, it follows from (4) that 

d j R=-R5 j p, (7) 

and it follows from this that djQ = —RSj Q. Summing over j, adding to (6) 
and evaluating at give 

T>Q(Zk) = P( Zk)-Q(Zk)Qu + [T, Q(Zk)]- 

2 Because of the fact pLX = R and our interpretation of Rij(x, £j) as Rij(x, wc 
are able to write R S p in place of pLS p. 

3 The meaning of 5 here and later is this: If U and V are matrix functions then U SV 
is the matrix with i,j entry ^ fc Uik(^k) Vkj(£k)- Thus RSQ is the matrix function with 
i,j entry ^2 k Rik{x, Cfc) Qkj (Cfc)- This makes it compatible with our use of S also as a 
multiplication operator so that, for example, (RS p) (A) = RS(pA). 
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If we define pij = Pij(£i) then we obtain 

Vq =p-q®u+ [t, q}. (8) 

Next we use the facts that D 2 — M commutes with L and that M com- 
mutes with X. It follows that 

[D 2 -M, K] = [D 2 -M, LX] =L[D 2 -M, X] =L[D 2 , X] =L(SD + D5). 

It follows from this that 

[D 2 -M, p] = pL5Dp + P LD5p. 

Applying both sides to A and using the fact that (D 2 — M)A = we obtain 

Q"(x) -xQ(x) = pLSQ' + pLDSQ. (9) 

The first term on the right equals R5Q'. For the second term observe that 

P LDX = pLXD + pL[D, X] = RD + pL5, 

so we can interpret that term as —R y 5Q (the subscript denotes partial 
derivative) where —R y (x, y) is interpreted as not containing the delta- function 
summand which arises from the jumps of R. With this interpretation of R y 
we can write the second term on the right as — R y 5Q. Thus, 

Q"{x) - x Q{x) = RSQ' - R y 5Q. 

Using this we obtain from (6) 

P' = xQ{x) + R5Q' - RySQ + Q'Ou- R X 5Q - [r, Q'], 

and then from (6) once more 

P' = x Q(x) + R 5 (P - Q 6 u + R 5 Q + [r, Q] ) - R y 5 Q 

+(P-Q@u + RSQ+[t, Q])Qu-R x 5Q-[t, P-QQu + R5Q + [t, Q]}. 
It follows from (5) that 

R x + Ry = -Q(x)GQ(y) + R5R+[t, p}. 

(We replaced RS p by RS R since, recall, R y does not contain delta-function 
summands.) We use this and also the identity R5[t, Q] — [r, RSQ] = — [r, R5]Q, 
and the fact that 5 and r commute. The result is that 
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P' = x Q(x) + R5P + Q{x) 6 Q5Q + (P - QOu+[t, Q])&u 

-2[r, P]<5Q- [r, P-Q6u+ [r, Q]]. 

It follows from (7) that djP = —R$j P- Summing over j, adding to the 
above and evaluating at £ k give 

V P(&) = ti k Q{ik) + Q(&) 6 + - Q(&) 9 u + [r, Q(&)]) 9 « 

-2[r, P(&, -^g-Ir, P(6)-Q(6)©«+[r, Q(&)]]. 
Hence Dp is equal to 

£g + qOqq + (p - q®u+ [r, g])6u - 2 [r, r] g - [r, p- gOn+ [r, g]]. 

Equivalently, in view of (8), 

Vp = £q + q&qq + Vq@u-2[T, r]q- [r, V q}. (10) 

Let us compute Vu. We have 
and so 

dkUij = -5 ik J Ai(£ k )p kj (£ k ,y)Ai(y)dy 

~ J J Ai(x) X»(a;) [R ik (x,£ k ) p kj (£ k ,y)] Ai(y) dx dy, 
where we use (7) again. This is equal to 

-<5 ifc Ai(£ fc ) Q kj {£ k ) - (QikiZk) ~ &ik Ai(6;)) Qkjitk), 

and so 

d k Uij = -Q ik {£k) Qkj{ik)- (11) 

This gives 

Vu=-qq. (12) 
Next, we find from (7) and (5) that 

VRfo, &) = -Q&)OQ(&) + [r, P(&, &)]■ 

This gives Pr = -g6g+[r, r], which is equation (3). 



To get equation (1) we apply V to (8) and use (10) and (12). We find 
that 

T> 2 q = ^q + q@ qq + T> q@u — 2 [r, r]q—[r, T> q\—T> qQu + qQ qq+[r, Vq] 

= Zq + 2qQqq-2[T, r] q, 

which is (1). 

Finally, to get equation (2) we use the fact that Xj(y) Pjk(y, x) is equal to 
Xk(x) times p' k -{x,y), where p' is the resolvent kernel for the matrix kernel 
with i,j entry Lji(x,y) Xj(y). Hence Qjk(x) is equal to Xk{x) times the 
Qkj(x) associated with Lji. Consequently for all the differentiation formulas 
we have for the Qkj(£,k), etc., there are analogous formulas for the Qjk(£,k), 
etc.. The difference is that we have to reverse subscripts and replace r by 
r* and r by — r. The upshot is that, by computations analogous to those 
used to derive (1), we derive another equation which can be obtained from 
(1) by making the replacements q — > q f , q — > q l , r — > r*, r — > — r and then 
taking transposes. The result is equation (2). 
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